Density fluctuations in the presence of spinodal instabilities 
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Density fluctuations resulting from spinodal decomposition in a non-equilibrium first-order chiral 
phase transition are explored. We show that such instabilities generate divergent fluctuations of 
conserved charges along the isothermal spinodal lines appearing in the coexistence region. Thus, 
divergent density fluctuations could be a signal not only for the critical end point but also for the 
first order phase transition expected in strongly interacting matter. We also compute the mean-field 
critical exponent at the spinodal lines. Our analysis is performed in the mean-field approximation 
to the NJL model formulated at finite temperature and density. However, our main conclusions are 
expected to be generic and model independent. 
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One of the central questions addressed in the con- 
text of QCD is the phase structure and the phase di- 
agram of strongly interacting matter at finite temper- 
ature and baryon number density [![. Based on cal- 
culations in effective models and on universality argu- 
ments one finds that the order of the transition from the 
hadronic to the quark gluon plasma phase depends on the 
number of quark flavors and on the value of the quark 
masses @, H 0, i, H 0, 1, S ■ For physical values of 
the parameters one expects that the transition at high 
temperature and low net baryon number density is con- 
tinuous. In the opposite limit of low temperature and 
large density the QCD phase transition is expected to be 
first order. This suggests that the phase diagram exhibits 
a critical end point (CEP), where the first order chiral 
transition of QCD terminates [1, EH- Recent, however 
still preliminary, results obtained in first principle calcu- 
lations of QCD, in lattice gauge theory (LGT), confirm 
the existence of such a point in the temperature T and 
chemical potential \x q plane [H, El, E3 ■ 

The search for the CEP has recently attracted con- 
siderable attention. It is of particular interest to iden- 
tify the position of the critical end point in the phase 
diagram and to study generic properties of thermody- 
namic quantities in its vicinity. The qualitative behavior 
of physical observables and their dependence on thermo- 
dynamic variables in the critical region can be studied in 
effective chiral models. However, a quantitative descrip- 
tion of the thermodynamics near the phase transition can 
only be obtained theoretically by solving QCD ab initio 
in LGT or phenomenologically in the context of experi- 
mental studies of heavy ion collisions. 

To locate the phase transition line in the QCD phase 
diagram one needs observables that are sensitive probes 
of the critical structure. Modifications in the magni- 
tude of fluctuations or the corresponding susceptibilities 
have been suggested as a possible signal for deconfinc- 
ment and chiral symmetry restoration [1, [ll|. In this 
context, fluctuations related to conserved charges are of 
particular interest [l5| . The fluctuations of baryon num- 



ber and electric charge diverge at the critical end point 
while they are finite along the cross over and first order 
phase boundaries [1, H, [la ] . Consequently, singular fluc- 
tuations of baryon number and electric charge as well as 
a non monotonic behavior of these fluctuations as func- 
tions of the collision energy in heavy ion collisions have 
been proposed as p ossible signals for the QCD critical 
end point [H, H, EH- However, the finiteness of the fluc- 
tuations along the first order transition depend on the 
assumption that this transition appears in equilibrium. 

A first order phase transition is intimately linked with 
the existence of a convex anomaly in the thermody- 
namic pressure 17|, which can be uncovered only in 
non-equilibrium systems. There is an interval of energy 
density or baryon number density where the derivative 
of the pressure is positive, dP/dV > 0. This anoma- 
lous behavior characterizes a region of instability in the 
(T, n 9 )-plane, where n q is the net quark number density. 
This region is bounded by the spinodal lines, where the 
pressure derivative with respect to volume vanishes. The 
derivative taken at constant temperature and that taken 
at constant entropy define the isothermal and isentropic 
spinodal lines, respectively. 

The consequences of spinodal decomposition have been 
discussed in connection with the chiral/ deconfinement 
phase transition in heavy ion collisions [l7], EH, El, HI, HH, 
|22| . Furthermore, spinodal decomposition plays a crucial 
role in the description of the 1st order nuclear liquid-gas 
transition in low energy nuclear collisions [13, [2l[ . It has 
also been argued that in the region of phase coexistence, 
a phase separation can lead to an enhancement of baryon 
[l8[ and strangeness fluctuations [h|. 

In this letter we consider the fluctuations of conserved 
charges along the spinodal lines, expected at finite net 
baryon density in the QCD phase diagram. We show that 
if the chiral phase transition is first order, then the fluctu- 
ations of the net densities of baryon number and electric 
charge diverge along the isothermal spinodal lines. Con- 
sequently, large fluctuations of these quantities may be 
a signal for a first order phase transition in the QCD 
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medium. We also compute the mean-field critical expo- 
nent of the quark susceptibility at the spinodal lines and 
compare with that at the CEP. 

In our study of fluctuations across the first order chi- 
ral phase transition we adopt the Nambu-Jona-Lasinio 
(NJL) model. For two quark flavors and three colors the 
NJL Lagrangian reads 



quark n; 



C = ip(i@ - m + n"io)i> + G s (^ip) 2 + (^if^ip) , (1) 

where m = diag(m u , mj) is the current quark mass, 
fx = diag(/i u , fid) are the quark chemical potentials and 
r are Pauli matrices. The strength of the interactions 
between the constituent quarks is controlled by the cou- 
pling G5A 2 = 2.44 with the three momentum cut-off 
A = 587.9 McV, introduced to regularize the ultraviolet 
divergences. The parameters are fixed to reproduce the 
pion decay constant in vacuum and the pion mass for 
m u = irid = 5.6 MeV. 

In the mean field approximation the thermodynamics 
of the NJL model is, for an isospin symmetric system, 
given by the thermodynamic potential [24j | : 



Q{T,fji;M)/V 



(M 
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d 3 p - 



E{p) 



AGs "J (2tt) 3 
- Tln(l - n<+>(£T,A0) -Tln(l - n^(f,T,n)] (2) 

where M = m — 2Gs(4>ip) is the dynamical quark 
mass, E(p) = y/p 2 + M 2 is the quasiparticle energy 

and (p, T,fx)= (l + exp [{E(p) =F fj) /T] ) 1 arc the 

quark/antiquark distribution functions. 

The dynamical mass M is obtained self consistently 
from the stationarity condition dCl/dM = 0, which im- 
plies 



M = m + 24G. 



d 3 p M 



(2tt) 3 E 



,(+) 
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(3) 



The relevant thermodynamic observables, the quark 
number density and the corresponding susceptibility, are 
given by 



dfi 
dfi ' 



dn q 
dfi 



(4) 



The NJL model has a generic, QCD like, phase dia- 
gram. It exhibits a critical end point that separates the 
cross over from the first order chiral phase transition. 
The relevant part of the phase diagram in the (T,n q )— 
plane is shown in Fig. 1. If the first order phase tran- 
sition takes place in equilibrium, there is a coexistence 
region, which ends at the critical end point. However, 
in a non-equilibrium first order phase transition, the sys- 
tem supercools/superheats and, if driven sufficiently far 
from equilibrium, it becomes unstable due to the convex 
anomaly in the thermodynamic pressure. In other words, 
in the coexistence region there is a range of densities and 
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FIG. 1: The phase diagram of the NJL model in the tem- 
perature (T, n q )-plane. The filled point indicates the CEP. 
The full lines starting at the CEP represent boundary of the 
coexistence region in equilibrium. The broken-curves are the 
isothermal whereas the dotted ones are the isentropic spinodal 
lines. 



temperatures, bounded by the spinodal lines, where the 
spatially uniform system is mechanically unstable. The 
location of the spinodal lines is determined by the condi- 
tions 



dP 
dV 



and 



zap 



0. 



(5) 



for the isothermal and isentropic spinodal lines, respec- 
tively. Both these lines are shown in Fig. 1. From the 
thermodynamic relation 



dP 
dV 



dP 
dV 



T 



dP 
df 



(6) 



it is clear that the isentropic spinodal lines are located 
inside the instability region bounded by the isothermal 
spinodals and that the two sets of lines coincide at T = 0. 
In the mean field approximation the isothermal instabil- 
ities disappear at the CEP, where the first order transi- 
tion ends, while the isentropic spinodal curves join below 
the CEP. We note that the isentropic spinodal lines will 
probably be modified considerably when fluctuations are 
properly included. This is because the specific heat Cy 
diverges at the CEP [l|, while in the mean-field approx- 
imation it remains finite. It then follows from (JB|) that 
both the isothermal and the isentropic pressure deriva- 
tives in ([5]) vanish at the CEP. Consequently, in a more 
complete description, the isentropic spinodal lines are ex- 
pected to move up in temperature and to also join at the 
CEP. 

Since the isothermal spinodal lines join at the CEP, as 
shown in Fig. 1, it is natural to explore how the charge 
fluctuations develop when going beyond the critical end 
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FIG. 2: The net quark number density fluctuations Xm/j/A 2 
as a function of the quark number density n q across the first 
order phase transition. The susceptibility Xmm was computed 
in the NJL model along the line of constant temperature T = 
50 MeV. 
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FIG. 3: The net quark number susceptibility Xmm/A in the 
stable and meta stable regions, as a function of the quark 
number density n q and temperature T. The full line corre- 
sponds to the critical temperature Tcep, and the dashed lines 
to the temperatures 30, 40, 50, 60, 70, 75 and 80 MeV. The 
dotted line shows the projection of the isothermal spinodal 
onto the (T,n q ) plane. 



point into the first order phase transition. In Fig. 2 we 
show the evolution of the net quark number fluctuations 
along a path of fixed T = 50 MeV in the (T, n q ) -plane. 
When entering the coexistence region, there is a singu- 
larity in Xfi.fi, that appears when crossing the isothermal 
spinodal lines, where the fluctuations diverge and the sus- 
ceptibility changes sign. Between the spinodal lines, the 
susceptibility is negative. This implies an instability of 
the baryon number fluctuations when crossing the tran- 
sition between the chirally symmetric and broken phases. 
An estimate of the growth rate for unstable density 



fluctuations, assuming a typical wave vector q = 100—300 
MeV and using the results of [25[ for the collisionless 
regime, yields a characteristic growth time r ~ 1 — 10 fm 
in agreement with [2(|. We note that the detailed calcu- 
lations of ref. [25| show that the growth rates of the most 
unstable mode deep inside the spinodal region depends 
only weakly on the collision rate. Thus, our estimate 
is expected to be valid also in the intermediate and hy- 
drodynamic regimes. Within the growth time an initial 
fluctuation grows by a factor e. Since this time is smaller 
than or on the order of the typical expansion time scale 
in relativistic heavy-ion collisions, we conclude that the 
spinodal instability may lead to observable fluctuations. 
A more quantitative estimate of the expected fluctuations 
in heavy-ion collisions would require a kinetic calculation, 
e.g. along the lines of ref. [29j . 

The behavior of Xp.p. seen in Fig. 2 is a direct conse- 
quence of the thermodynamic relations 



dP 
W 

dP 
dV 



n 2 q 1 



n 2 q Xtt 



'-X/J.T 



Xfifj. 



V 



XfJ.fJ.XTT ~ XfiT 



(7) 
(8) 



which connect the pressure derivatives with the suscep- 
tibilities Xxy = —d 2 Q,/dxdy. Along the isothermal spin- 
odal lines the pressure derivative in |(7j) vanishes. Thus, 
for non- vanishing density n q , Xfifi must diverge to satisfy 
(J7]). Furthermore, since the pressure derivative dP/dV\T 
changes sign when crossing the spinodal line, there must 
be a corresponding sign change in Xfifn as seen in Fig. 2. 
Due to the linear relation between Xmmi the isovector sus- 
ceptibility xi an d the charge susceptibility XQ H3i the 
charge fluctuations are also divergent at the isothermal 
spinodal line. Thus, in heavy-ion collisions, fluctuations 
of the baryon number and electric charge could show en- 
hanced fluctuations, as a signal of the spinodal decom- 
position. The spinodal phase separation can also lead to 
fluctuations in strangeness [l9[ and isospin densities [3(| • 

At the isentropic spinodal line the baryon number sus- 
ceptibility is in general finite. This is also true for the 
other susceptibilities appearing in Eq. ([5]). The isen- 
tropic spinodal line ([S]) corresponds to a zero of the nu- 
merator in {8j and of the velocity of hydrodynamic sound 
waves. In the hydrodynamic limit the instability sets in 
at the isothermal spinodal, but the growth rate becomes 
large only at the isentropic spinodal line [25[ . 

In the case of an equilibrium first order phase tran- 
sition, the density fluctuations do not diverge; in the 
coexistence region the susceptibility is a linear combi- 
nation of the positive susceptibilities above and below 
the phase boundary. Thus, the fluctuations increase as 
one approaches the CEP along the first order transition 
and decrease again in the cross over region. This led 
to the prediction of a non-monotonous behavior of the 
fluctuations with increa sing beam energy as a signal for 
the existence of a CEP [llf, [28| . We stress that strictly 
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speaking this is relevant only for the idealized situation 
where the first order phase transition takes place in equi- 
librium. In the more realistic non-equilibrium system one 
expects fluctuations in a larger region of the phase dia- 
gram, i.e., over a broader range of beam energies, due to 
the spinodal instabilities. 

The critical exponent at the isothermal spinodal line 
is found to be *y = 1/2, with X/x/x ~ 

(fi - n c ) 7 , while 

7 = 2/3 at the CEP, in agreement with the mean-field 
results Thus, the singularities at the two spinodal 

lines conspire to yield a somewhat stronger divergence as 
they join at the CEP. The exponents are renormalized by 
fluctuations, but the smooth evolution of the singularity 
from the spinodal lines to the CEP, illustrated in Fig. 3, 
is expected to be generic. 

We have shown that the net quark number fluctuations 
diverge at the isotherm al sp inodal lines of the first order 
chiral phase transition 3j|. As the system crosses this 
line, it becomes unstable with respect to spinodal de- 
composition. The unstable region is in principle reach- 
able in non-equilibrium systems, created e.g. in heavy 
ion collisions. This means that large fluctuations of the 
density is expected not only at the second order critical 
end point but also at a non-equilibrium first order tran- 
sition. In fact, the signal from the first-order transition 
may be much stronger than that from the CEP. Model 
calculations suggest that the critical region of enhanced 
susceptibility around the CEP is fairly small [§, [lj| [28[ , 



while it is large in the spinodal region, where the fluctu- 
ations appear due to the divergence of and due to 
the mechanical instability of the system. We stress that 
there is a close relation between the singularities at the 
CEP and the spinodal lines. 

The properties of different susceptibilities were ob- 
tained within the NJL model in the mean field approxi- 
mation. However, the singular properties of charge sus- 
ceptibilities in the presence of spinodal instabilities are 
quite general. They appear due to the straightforward 
thermodynamic relation between the pressure derivatives 
and charge susceptibilities. Thus, although the NJL 
model is non-confining, the results presented in this pa- 
per are expected to be robust on a qualitative level. 
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